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Abstract. The ANTARES telescope is well-suited for detecting astrophysical transient neu¬ 
trino sources as it can observe a full hemisphere of the sky at all times with a high duty cycle. 
The background due to atmospheric particles can be drastically reduced, and the point-source 
sensitivity improved, by selecting a narrow time window around possible neutrino produc¬ 
tion periods. Blazars, being radio-loud active galactic nuclei with their jets pointing almost 
directly towards the observer, are particularly attractive potential neutrino point sources, 
since they are among the most likely sources of the very high-energy cosmic rays. Neutrinos 
and gamma rays may be produced in hadronic interactions with the surrounding medium. 
Moreover, blazars generally show high time variability in their light curves at different wave¬ 
lengths and on various time scales. This paper presents a time-dependent analysis applied 
to a selection of flaring gamma-ray blazars observed by the FERMI/LAT experiment and 
by TeV Cherenkov telescopes using five years of ANTARES data taken from 2008 to 2012. 
The results are compatible with fluctuations of the background. Upper limits on the neu¬ 
trino fluence have been produced and compared to the measured gamma-ray spectral energy 
distribution. 
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1 Introduction 

Neutrinos are unique messengers for studying the high-energy Universe as they are neutral, 
stable, interact weakly, and travel directly from their sources without absorption or deflection. 
Therefore, the reconstruction of the arrival directions of cosmic neutrinos would allow both the 
sources of the cosmic rays - supernova remnant shocks, active galactic nuclei jets, gamma-ray 
bursts, etc. [ 1 ] - and the relevant acceleration mechanisms acting within them to be identified. 

The high-energy extragalactic sky is dominated by active galactic nuclei (AGN). Their 
spectral energy distribution can be described by two components: a low-energy one from 
radio to X-rays and a high-energy one from X-rays to very high-energy gamma rays. The 
low-energy component is generally attributed to synchrotron radiation in the relativistic jet 
by a non-thermal population of accelerated electrons and positrons; the origin of the second 
component is still under discussion. In leptonic models [2], it is ascribed to an inverse Compton 
process between the electrons and a low-energy photon field (their own synchrotron radiation, 
or external photons), while in hadronic models it originates from synchrotron emission by 
protons and secondary particles coming from p -7 or p-p interactions [3, 4]. Associated with 
these very high-energy gamma rays from decays, the decay of the charged pions gives rise 
to a correlated neutrino emission. 

Flat-spectrum radio quasars (FSRQs) and BL Lacs, together classified as blazars, exhibit 
relativistic jets pointing almost directly towards the Earth, and are among the most violent 
variable high-energy phenomena in the Universe [5]. Blazars are particularly attractive po¬ 
tential neutrino sources, since they are among the most likely sources of the very high-energy 
cosmic rays [ 6 , 7]. Several theoretical models predict high-energy neutrino emission from 
blazars that yield different shapes and normalisations for the expected energy spectrum [ 8 - 
11]. In addition, the models also suffer from large uncertainties that originate, for instance, 
from unknowns in the model parameters, the luminosity functions and the source evolution. 
For example, the FSRQs are predicted to be more promising neutrino candidates than BL 
Lacs in Ref. [12], while in Ref. [13] the opposite is predicted. Several authors even estimate 
more optimistic spectra with spectral indexes up to one [14, 15]. The E~‘^ spectrum, generally 
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expected from Fermi acceleration of cosmic rays in astrophysical sources, is used as a reference 
spectrum. An energy cutoff seems to be present in most sources observed in gamma rays [16]. 
For these reasons, to cover the majority of the range allowed by the models accessible to the 
ANTARES sensitivity (see Section 2), four neutrino-energy spectra are tested in this analysis: 
E~'^, E~‘^ exp{—E/lO TeV), E~‘^ exp{—E/I TeV) and E~^, where E is the neutrino energy. 

In the ANTARES telescope [17], events are primarily detected underwater by observing 
the Cherenkov light induced by relativistic muons in the darkness of the deep sea. Owing to 
their low interaction probablility, only neutrinos have the ability to cross the Earth. Therefore, 
an upgoing muon is an unambiguous signature of a neutrino interaction close to the detec¬ 
tor. The detection of cosmic neutrinos with neutrino telescopes is very challenging because 
of the small neutrino interaction cross-section and the high background of atmospheric neu¬ 
trinos from cosmic-ray interactions in the atmosphere. To distinguish astrophysical neutrino 
events from background events (muons and neutrinos) generated in the atmosphere, energy 
and direction reconstructions have been used in several searches [18-20]. To improve the 
signal-to-noise discrimination, the arrival time information can be used, significantly reduc¬ 
ing the effective background [21, 22]. Blazars are known to show time variability at different 
wavelengths and on various time scales [23-25]. The associated neutrino emission may ex¬ 
hibit similar variability, and this is used in time-dependent methods to improve the detection 
probability with respect to time-integrated approaches. 

In this paper, the results of a time-dependent search for cosmic neutrino sources using 
the ANTARES data taken from 2008 to 2012 is presented. This extends a previous ANTARES 
analysis [26] where only the last four months of 2008 were considered. The analysis is applied 
to a list of promising blazar candidates detected in GeV gamma rays by the LAT instrument 
onboard the FERMI satellite, and to a list of blazar flares reported by TeV gamma-ray 
experiments (H.E.S.S., MAGIC and VERITAS). After a brief description of the apparatus, 
the data selection and the detector performances are presented in Section 2. The point-source 
search algorithm used in this time-dependent analysis is explained in Section 3. The results of 
the GeV and TeV flare searches are presented in Sections 4 and 5, respectively, and discussed 
in Section 6. 

2 The ANTARES neutrino telescope 

The ANTARES collaboration completed the construction of a neutrino telescope in the 
Mediterranean Sea with the connection of its twelfth detector line in May 2008 [17]. The 
telescope is located 40 km off the Southern coast of France (42°48’N, 6°10’E), at a depth of 
2475 m. It comprises a three-dimensional array of 10" photomultipliers [27], each housed in 
a glass sphere (Optical Modules, OMs [28]), distributed along twelve slender lines anchored 
at the sea bottom and kept taut by a buoy at the top. The lines are connected to a central 
junction box, which in turn is connected to shore via an electro-optical cable. Since lines 
are subject to the sea current and can change shape and orientation, a positioning system 
comprising hydrophones, compasses and tiltmeters is used to monitor the detector geome¬ 
try [29] . The main goal of the experiment is to search for neutrinos of astrophysical origin by 
detecting high-energy muons (>100 GeV) created by neutrino charged-current interactions in 
the vicinity of the detector. 

The arrival time and intensity of the Gherenkov light on the OMs are digitised into 
‘hits’ and transmitted to shore, where events containing muons are separated from the optical 
backgrounds due to natural radioactive decays and bioluminescence, and stored on disk. A 
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detailed description of the detector and of the data acqnisition is given in Ref. [17, 30]. The 
data nsed in this analysis were taken with the fnll detector dnring the period from September 
6th, 2008 np to December 31st, 2012 (54720-56292 modified Jnlian day). Filters are applied 
in order to exclnde periods in which the biolnminescence-indnced optical backgronnd was 
high. The resnlting effective livetime is 1044 days. 

Atmospheric nentrinos are the main sonrce of backgronnd in the search for astrophysical 
nentrinos. These nentrinos are prodnced from the interaction of cosmic rays in the Earth’s 
atmosphere. To acconnt for this backgronnd, mnon nentrino events are simnlated nsing the 
GENHEN package [31] according to the parametrisation of the atmospheric nentrino flnx 
from Ref. [32] . An additional sonrce of backgronnd is dne to mis-reconstrncted atmospheric 
mnons. Downgoing atmospheric mnons are simnlated with the program MURAGE [33, 34] 
which provides mnon bnndles at the detector. The fnll Monte Garlo chain, which inclndes 
the simnlation of Gherenkov photon prodnction and of the detector response, is described in 
Ref. [35]. 

The applied track reconstruction algorithm is the same as that nsed in the standard 
ANTARES point-sonrce search [19]. This algorithm derives the mnon track parameters that 
maximise a likelihood fnnction bnilt from the difference between the expected and the mea- 
snred arrival times of the hits from the Gherenkov photons emitted along the mnon track. 
This maximisation takes into acconnt the Gherenkov photons that scatter in the water and the 
additional photons that are generated by secondary particles (e.g. electromagnetic showers 
created along the mnon trajectory). 

The mnon track reconstrnction retnrns two qnality parameters, namely the track-ht 
qnality parameter. A, and the estimated angnlar nncertainty on the htted mnon track direc¬ 
tion, j3. Gnts on these parameters are nsed to improve the signal-to-noise ratio. To ensnre a 
good directional reconstrnction of the selected nentrino candidates, /? < 1° is reqnired. Fig- 
nre 1 shows the distribntion of the error estimate (3. This cnt rejects 47% of the atmospheric 
mnons which are mis-reconstrncted as npgoing tracks. After qnality cnts, all events with a 
zenith angle cos(0) > -0.15 are selected. Gompared to a strict npgoing event selection, this 
cnt provides a gain in visibility of on average 15% for all sonrces with a declination above 
-20°. Fignre 2 shows the distribntion of the reconstrncted cosine of the zenith angle for both 
data and simnlation. The valne of the cnt on A is optimised for each sonrce on the basis 
of maximising a model discovery potential [22] for a 3a signihcance level for each nentrino 
spectrnm. The distribntion of A for events selected after applying the /3 and cos(0) cnts is 
shown in Fignre 3. The optimnm A valnes range from -5.5 to -5.0 depending on the sonrce 
and the backgronnd characteristics dnring the flares. 

The energy of each events is estimated by exploiting the correlation between the energy 
deposition, dE/dX, and the primary energy [36, 37]. The systematic nncertainty in the dE/dX 
energy estimator is 10%. This is acconnted for by smearing the simnlated signal valnes by a 
ganssian fnnction with this RMS valne. 

The pointing accnracy has been determined directly from the data nsing the moon 
shadow [38] and is of the order of 0.63° (G.L. 90%). Fignre 4 shows the cnmnlative distribntion 
of the angnlar difference between the reconstrncted mnon direction and the nentrino direction 
with an assnmed spectrnm proportional to E~‘^, E~‘^ exp(—Fl/10 TeV), E~‘^ exp(—Fl/1 TeV) 
and E~^. For example, the median resolntion is estimated to be 0.43±0.1° for the reference 
E~‘^ energy spectrnm. 

The acceptance (with nnits of GeV cm^ s^) is the proportionality factor between a 
given energy flnx and the corresponding nnmber of signal events expected in the detector 
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Figure 1. Distribution of the estimated error on the direction of the reconstructed muon track after 
applying a zenith angle cut cos(6*) > -0.15 and a cut on the quality variable A > -5.3. The blue line 
shows the simulated upgoing atmospheric neutrinos, the red line the mis-reconstructed atmospheric 
muons, the green line the sum of both contributions, and the black crosses the data. The vertical 
dashed line with the arrow shows where the selection cut is applied (/3 < 1°). 


after selection. Thus, a certain neutrino spectrum has to be assumed. For example, for 
A > —5.3, the acceptances as a function of the source declination are shown for each considered 
spectra in Figure 5. For the limit setting, a 15% systematic uncertainty on the acceptance is 
considered [19]. 

3 Time-dependent search method 

The time-dependent point-source analysis is performed using an unbinned method based on 
a likelihood-ratio maximisation. The data are parameterised as a two-component mixture of 
signal and background. The signal is expected to be small so that the full data sample (N 
events) can be treated as background. The expected number of events are A/s (unknown) and 
A/b (known) and the probabilities for signal and background for an event i, at time ti, energy 
estimate dE/dXi, declination di and 'I'j are Si and Bi respectively. S and B respectively. The 
probability V and the likelihood L are: 


'Pi = A/s5i -b M-qBi 


(3.1) 


ln£ 


N 


M^Bi 


[A/s -b A/b] 


(3.2) 


\ 2=1 / 

To discriminate the signal-like events from the background ones, these probabilities are 
described by the product of three components related to the direction, energy, and timing of 
each event. For an event z, the signal probability is: 
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Figure 2. Distribution of the reconstructed cosine of the zenith angle of the events (black crosses) 
with /3 < 1° and A > -5.3. The simulated distributions are shown for atmospheric muons (red) and 
upgoing neutrinos (blue) and the green line is the sum of both components. The vertical dashed line 
with the arrow shows where the cut on the zenith angle is applied in order to select mainly upgoing 
events (those with larger zenith angles). 


^ . s---’^sy{dE/dXi) ■ + lag) (3.3) 

where is a parameterisation of the point spread function, i.e., 5®P^“('E'i(Q;s, <5s)) the 

probability to reconstruct an event i at an angular distance from the true source location 
(q!s,(5s). The energy PDF is parametrised with the normalised distribution of the muon 

energy estimator of an event according to the studied energy spectrum. Figure 6 illustrates 
the energy PDF of these four tested energy spectra. For the generation of signal events, the 
correlation between the angular error and the energy of the muon track is taken into account. 
The shape of the time PDF, 5*™®, for the signal event is extracted directly from the gamma- 
ray light curve assuming the proportionality between the gamma-ray and the neutrino fluxes. 
A possible lag of up to ±5 days has been introduced in the likelihood to allow for small lags 
in the proportionality. This corresponds to a possible shift of the entire time PDF. The lag 
parameter is fitted in the likelihood maximisation together with the number of fitted signal 
events in the data. 

The background probability for an event i is: 

Bi = i3®P'^®"((5i) • B^^^^^y{dE/dXi) • i3'“"(ti) (3.4) 

The directional PDF the energy PDF and the time PDF for the 

background are derived from data using, respectively, the observed declination distribution of 
selected events in the sample, the measured distribution of the energy estimator, and the ob¬ 
served time distribution of all the reconstructed muons. Figure 7 shows the time distribution 
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Figure 3. Distribution of the reconstruction quality variable A for tracks (black crosses) which have 
an angular uncertainty estimate ^ < 1° and zenith angle cos(6*) > -0.15. The simulated distributions 
are shown for atmospheric muons (red) and upgoing neutrinos (blue) and the green line is the sum of 
both components. The vertical dashed area with the arrow shows where the optimised selection cuts 
stand for the various tested sources. 


of all reconstructed events, including both downgoing and selected events. Once normalised, 
the hrst downgoing distribution is used directly as the time PDF for the background. Null 
values indicate the absence of data taken during these periods (e.g. detector in maintenance) 
or data with a very poor quality (high bioluminescence). 

The goal of the unbinned search is to determine, in a given direction in the sky and at 
a given time, the relative contribution of each component, and to calculate the probability to 
have a signal above a given background model. This is done via the test statistic. A, dehned 
as the ratio of the probability for the hypothesis of background and signal (Hgig+hkg) over the 
probability of only background (Ffbkg): 


, _ I 'Pjxil P^sig-l-bkg (A/s ) ) 

- ^ ^ P(x.lifbkg) 


(3.5) 


where A/s and N are respectively the unknown number of signal events and the total number 
of events in the considered data sample, and Xi are the observed event properties (hj, RAi, 
dE/dXi and tj) 

The evaluation of the test statistic is performed by generating pseudo-experiments sim¬ 
ulating background and signal in a 30° cone around the considered source according to the 
background-only and background plus signal hypotheses. The direction of the background 
events are generated by randomly sampling the declination distribution from the background 
probablility B and the right ascension from a uniform distribution. Signal events are simulated 
by hrst sampling the time and the energy of the simulated events by randomly generating 
from and respectively. Then, the angular distance from the coordinates of the 
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Figure 4. Cumulative distribution of the angle, a, between the true Monte Carlo neutrino direction 
and the reconstructed muon direction for neutrino events selected for this analysis with E~^ (black 
line, with median 0.43°), E~'^ exp{—E/10 TeV) (red line, with median 0.55°), E~'^ exp(—E/l TeV) 
(green line, with median 0.90°) and E~^ (blue line, with median 0.29°). 


studied source is sampled as a function of the declination and of the estimated 

energy. 

The null hypothesis is given by A/s = 0 (Aq) (i.e. the background-only hypothesis i/fbkg)- 
The obtained value of A for the data, Adataj is then compared to the distribution of Aq obtained 
by pseudo-experiments. Large values of Adata compared to the distribution of Aq reject the 
null hypothesis with a conhdence level depending on the fraction of the Aq distribution above 
^data- This fraction of trials above Adata is referred to as the p-value. The discovery potential 
is then dehned as the average number of signal events required to achieve a p-value lower 
than 2.7 10“^ (5.7 10“'^) (3(5)cr) in 50% of the trials. 

The performance of the time-dependent analysis is computed with a toy experiment with 
a source assuming a square-shaped flare with a width varying from 1 to 2000 days assuming 
a flat background period of 2000 days. Figure 8 shows the average number of events required 
for a 5cJ discovery for a single source located at a declination of -40° and following an E~‘^ 
energy spectrum, as a function of the total width of the flare periods. These numbers are 
compared to those obtained without the selection of time intervals corresponding to flaring 
periods. For time ranges characteristic of flaring activity, the time-dependent search presented 
here improves the discovery potential by on-average a factor 2-3 with respect to a standard 
time-integrated point-source search [19] under the assumption that the neutrino emission is 
correlated with the gamma-ray flaring activity. 
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sin(5) 

Figure 5. Acceptance of the analysis as a function of the sine of the declination for E~'^ (black), 
E~‘^ exp{—E/10 TeV) (red), E~‘^ exp{—E/l TeV) (green) and E~^ (blue) neutrino energy spectra. 
The events are selected with /? < 1°, cos(0) > -0.15 and A > -5.3. The different shape between 
the acceptance for E~^ and the rest of the neutrino spectra is due to the Earth absorption for very 
high-energy events. 


4 Search for neutrino emission from gamma-ray flares detected by FERMI 

The time-dependent analysis described in the previous section is applied to bright and variable 
Fermi blazar sources reported in the second Fermi LAT catalogue [39] and in the LBAS 
catalogue (LAT Bright AGN sample [16]). The sources located in the part of the sky visible 
to ANTARES ((5 < 35°) with a flux greater than 10“® photons • cm“^ • s“^ above 1 GeV, a 
test statistic TS > 25 (corresponding to a detection significance of more than 4 sigma) and 
a significant time variability are selected. This list is completed by adding sources reported 
as flaring in the Fermi Flare Advocates in 2011 and 2012 [40]. The hnal list includes a total 
of 153 sources. 

Light curves for the selected sources are produced using the Fermi Public Release Pass 
7 data using the source class event selection (evclass=2) and the Fermi Science Tools v9r35pl 
package [41]. Following the standard event selection cuts proposed by the Fermi-LAT Gollab- 
oration [42], the data are hltered using the gtselect tool to select only events which are most 
likely gamma-rays. Light curves are computed from the photon counting data in a cone of 
two-degree radius around each source direction, corrected by the total exposure. With this 
method the diffuse background contributions are not subtracted. This limits the validity of 
this method to only bright gamma-ray sources, but does not affect the flare identihcation 
in these sources. Therefore, sources which are close to the galactic plane (galactic latitude 
< 10°) or have others sources within a two-degree cone (or 3° for very bright sources) are 
excluded. The baseline is then removed for the time PDF dehnition. The resulting light 
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Figure 6. Distribution of the energy PDF of the background events (black) and of the signal events as¬ 
suming an energy spectrum following E~'^ (red), E~'^ exp{—E/lO TeV) (green), E~'^ exp{—E/I TeV) 
(blue) and E~^ (purple). Each distribution is normalised to unity. 


curves correspond to the one-day-binned time evolution of the average gamma-ray flux above 
a threshold of 100 MeV from August 2008 to December 2012. This method has the main 
advantage of producing continuous and complete gamma-ray light curves. Figure 9 (a) shows 
the resulting light curve for the source 3C273. This aperture photometry method agrees with 
the results of the likelihood analysis (an alternative method used for the LAT data analy¬ 
sis). The exposure is calculated using the gtexposure tool, which is also part of the Fermi 
framework. 

A maximum likelihood block (MLB) algorithm [43-45] is used to remove noise from the 
light curves by iterating over the data points and selecting periods during which data are 
consistent with a constant flux within statistical errors. The description of the light curve in 
terms oi Nb periods of constant emission is obtained with the following likelihood: 


InL = 



i&riB 


(4.1) 


where Xi is the flux measurement of the data point i and ai is the error of the flux 
measurement i. Since this likelihood would maximise when there are as many periods as data 
points, while only a few periods are required to describe the light curve, the following prior 
V is added to the likelihood: 

\nV = -NslnAf 
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Figure 7. Time distribution of the reconstructed events. Upper histogram (black): distribution of 
well-reconstructed events (including downgoing muons). Bottom histogram (red): distribution of the 
events selected by this analysis. 



Figure 8. Average number of events required for a 5(t discovery (50% probability) for a source located 
at a declination of -40° and following an E~‘^ energy spectrum as a function of the total width of the 
flaring periods (solid line). These numbers are compared to those obtained without using the timing 
information (dashed line). 
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where J\f is the number of data points in the light curve. The algorithm searches in each 
iteration for the best two flux changes in the emission in each period already found and keeps 
them for the next iteration. The Figure 9 (b) shows the optimisation of the likelihood in the 
case of 3C273. 

High-state periods are dehned using a simple and robust method. The value of the 
steady state (i.e. baseline, BL) and its fluctuation {(Jbl) are determined with a Gaussian ht 
of the lower part of the distribution of the flux data points (Figure 9 (c)). The flaring periods 
are dehned in three main steps. Firstly, seeds are identihed by searching for points with an 
amplitude, or blocks with a huence above BL 5aBL- Then, each period is extended forward 
and backward up to an emission compatible with BL -|- laBL- An additional delay of 0.5 
days is added before and after the hare in order to take into account that the precise time 
of the hare is not known (one-day binned light curve). Finally, spurious hares are discarded 
by requiring that each hare is visible in the three light curves produced with a gamma-ray 
threshold of 100, 300 and 1000 MeV. Figure 9 (d) shows an example of one hare of 3C273. 
With the above dehnition, a hare has a width of at least two days. Under the hypothesis that 
the neutrino emission follows the gamma-ray emission, the signal time PDF is the normalised 
de-noised light curve with only the high state periods (the other periods are set to zero). The 
hnal list includes 41 bright and variable Fermi blazars: 33 FSRQs, 7 BL Lacs and 1 unknown 
identihcation. Figure 10 shows the position of the selected Fermi Blazars together with the 
ANTARES visibility. The main characteristics of these blazars are reported in Table 1. 

The results of the search is summarised in Table 2. Only three sources, 3C279, PKS10235- 
618 and PKS1124-186, have a pre-trial p-value lower than 10%. The lowest p-value, 3.3%, is 
obtained for 3C279 where one event is coincident with a large gamma-ray flare detected by 
Fermi/LAT in November 2008. Figure 11 shows the Fermi gamma-ray light curve of 3C279 
with the time of the neutrino events, the estimated energy distribution, and the angular 
distribution of the events around the position of this source. This coincident event is recon¬ 
structed with 89 hits (energy deposition dE/dX = 244 in arbitrary units) distributed on ten 
lines with a track ht quality A = —4.4. The particle track direction is reconstructed at 0.3° 
from the source location. This event has already been reported in the previous analysis [26]. 
The post-trial probability, computed by taking into account the 41 searches, is 67%, and is 
thus compatible with background huctuations. 

In the absence of a discovery, upper limits on the neutrino huence, at 90% conhdence 
level are computed using 5-95% of the energy range as: 


Ty = At 



dE-E 


diV 

dE 


The emission duration. At, is computed using the time PDF and the effective livetime. 
The limits include the systematic errors and are calculated according to the classical (frequen- 
tist) method for upper limits [46]. Figure 12 gives these upper limits. IceCube has performed 
a similar time-dependent analysis [47] using data from 2008 to 2012, with 19 sources in com¬ 
mon with the list presented in Table 1. For sources in the Southern Hemisphere, the limits 
computed by IceCube are of the same order of magnitude as the ANTARES limits while they 
are a factor 10 stronger for the sources in the Northen hemisphere thanks to IceCube’s larger 
detector volume. 
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Table 1. List of bright variable Fermi blazars with significant flares selected for this analysis. 


Name 

OFGL name 

Glass 

RA [°] 

Dec [°] 

Redshift 

3C 454.3 

J2254.0+1609 

FSRQ 

343.50 

16.15 

0.859 

4C +21.35 

J1224.9+2122 

FSRQ 

186.23 

21.38 

0.434 

PKS 1510-08 

J1512.7-0905 

FSRQ 

228.18 

-9.09 

0.360 

3C 279 

J1256.1-0548 

FSRQ 

194.03 

-5.8 

0.536 

PKS 1502+106 

J1504.3+1029 

FSRQ 

226.10 

10.49 

1.839 

PKS 2326-502 

J2329.2-4956 

FSRQ 

352.32 

-49.94 

0.518 

3C 273 

J1229.1+0202 

FSRQ 

187.28 

2.05 

0.158 

AO 0235+164 

J0238.6+1636 

BLLac 

39.65 

16.61 

0.940 

PKS 0426-380 

J0428.6-3756 

FSRQ 

67.17 

-37.93 

1.110 

4C +28.07 

J0334.3-3728 

FSRQ 

53.58 

-37.47 

1.206 

PKS 0454-234 

J0457.1-2325 

FSRQ 

74.28 

-23.43 

1.003 

PKS 1329-049 

J1332.0-0508 

FSRQ 

203.01 

-5.14 

2.150 

PKS 0537-441 

J0538.8-4405 

BLLac 

84.71 

-44.08 

0.896 

4C +14.23 

J0725.3+1426 

FSRQ 

111.33 

14.44 

1.038 

PMNJ 0531-4827 

J0532.0-4826 

UNID 

83.01 

-48.44 

/ 

PKS 0402-362 

J0403.9-3604 

BLLac 

60.99 

-36.07 

1.417 

PKS 1124-186 

J1126.6-1856 

FSRQ 

171.66 

-18.95 

1.048 

Ton 599 

J1159.5+2914 

FSRQ 

179.88 

29.25 

0.725 

PKS 2142-75 

J2147.4-7534 

FSRQ 

326.87 

-75.58 

1.139 

PKS 0208-512 

J0210.8-5100 

FSRQ 

32.70 

-51.2 

1.003 

PKS 0235-618 

J0237.1-6136 

FSRQ 

39.29 

-61.62 

0.467 

PKS 1830-211 

J1833.6-2104 

FSRQ 

278.41 

-21.08 

2,507 

PKS 2023-07 

J2025.6-0736 

FSRQ 

306.42 

-7.61 

1.388 

PKSB 1424-418 

LJ1428.0-4206 

FSRQ 

217.01 

-42.10 

1.522 

PMNJ 2345-1555 

J2345.0-1553 

FSRQ 

356.27 

-15.89 

0.621 

OJ 287 

J0855.4+2009 

BLLac 

133.85 

20.09 

0.306 

PKS 0440-00 

J0442.7-0017 

FSRQ 

70.69 

-0.29 

0.845 

PKS 0250-225 

J0252.7-2218 

FSRQ 

43.20 

-22.31 

1.419 

B22308+34 

J2311.0+3425 

FSRQ 

347.77 

34.43 

1.187 

B21520+31 

J1522.1+3144 

FSRQ 

230.54 

31.74 

1.487 

PKS 1730-13 

J1733.1-1307 

FSRQ 

263.28 

-13.13 

0.902 

PKS0244-470 

J0245.9-4652 

FSRQ 

41.06 

-47.06 

1.385 

PKS 0301-243 

J0303.4-2407 

BLLac 

45.87 

-24.13 

0.26 

CTA 102 

J2232.4+1143 

FSRQ 

338.12 

11.72 

1.037 

OG 050 

J0532.7+0733 

FSRQ 

83.19 

7.56 

1.254 

PMNJ 2331-2148 

J2330.9-2144 

FSRQ 

352.75 

-21.74 

0.563 

PKS0805-07 

J0808-0751 

BLLac 

122.06 

-7.85 

1.837 

PKS 2320-035 

J2323.6-0316 

FSRQ 

350.91 

-3.28 

1.411 

PKS 2227-08 

J2229.7-0832 

FSRQ 

337.44 

-8.55 

1.560 

OJ248 

J0830.5+2407 

FSRQ 

127.72 

24.18 

0.942 

PKS 2233-148 

J2236.5-1431 

BLLac 

339.13 

-14.53 

0.325 




Table 2. Results of the search for neutrinos in coincidence with Fermi blazar flares (three first rows) 
and TeV flares (last row). The total duration of all identified flares At, the optimised Aopt cuts, 
the number of required events for a 3a discovery (TVaa ), the number of fitted signal events by the 
likelihood (Nfn), the fitted time lag (Lag) and the corresponding pre-trial (post-trial) probability are 
given together with the energy spectra. 


Source 

At 

^opt 

Vs. 

Nf^t 

Lag 

P-value 

Post-trial 

Spectrum 

3G279 

279 d 

-5.3 

2.5 

0.8 

-4 d 

0.033 

0.67 

E-'^ 

PKS1124-186 

73 d 

-5.4 

3.1 

0.7 

-t-4 d 

0.059 

0.94 

E;-^exp(-E/lTeV) 

PKS0235-618 

25 d 

-5.7 

1.5 

0.6 

-4 d 

0.045 

0.91 

V-^exp(-E/10TeV) 

PKS0447-439 

10 d 

-5.4 

0.75 

0.1 

-1-5 d 

0.10 

0.55 

V-^exp(-E/lTeV) 


5 Search for neutrino emission from gamma-ray flares detected by TeV 
telescopes 

Ground-based observatories such as H.E.S.S., MAGIC and VERITAS cannot monitor sources 
continuously, because they generally have a reduced field of view and a low duty cycle (e.g. 
only moonless nights). Nevertheless, these telescopes detect photons with energies from a 
few hundred GeV to a few TeV that may be better correlated with the neutrinos to which 
ANTARES is sensitive. These observatories often emit alerts reporting flares to Astronomer’s 
Telegram or directly in a dedicated paper. When the start and stop times of the flare are 
known, the time PDF is assumed to be a single square-shaped flare with a minimum width 
of one day. Often, the beginning and the end of the flaring activity cannot be constrained 
accurately. In this case, a simple time cut is used, taking a time window including two 
days before and after the identified flare. Table 3 presents the list of seven TeV flares, 
their characteristics, and the publications from where this information is extracted. The 
flares are chosen for this analysis according to the same visibility criteria as for Fermi/LAT 
observations. The same analysis as described previously is performed assuming the same four 
energy spectra. 

Six of the seven flares tested show no excess of events in the vicinity of the corresponding 
sources in the selected time windows. Only the blazar PKS0447-439 shows a pre-trial p- 
value lower than 10% in the case of the assumed E~‘^ ex.p{—E/l TeV) energy spectrum (cf 
Table 2). The corresponding post trial p-value is 55%, and is also consistent with background 
fluctuations. In the absence of a signal, upper limits on the neutrino fluence at 90% confidence 
level are computed including the systematic errors (Figure 13). 

Table 3. List of TeV flares reported by the three telescopes H.E.S.S., MAGIC and VERITAS in the 
2008-2012 period. 


Name 

Telescope 

RA 

Dec 

Flaring days (MJD) 

Reference 

4G-f21.35 

MAGIG 

186.2 

21.4 

55364-5 

arXiv: 1101.4645 

PG 1553+113 

MAGIG 

239.0 

11.2 

55980-91 56037-8 

arXiv: 1109.5860 

PKS 1424+240 

MAGIG 

216.8 

23.8 

54940-60 

arXiv: 1109.5860 

lES 1218+30.4 

VERITAS 

185.4 

30.2 

54860-5 

arXiv: 1005.3747 

lES 0229+200 

VERITAS 

38.2 

20.3 

55118-31 

arXiv: 1307.8091 

J1943+213 

H.E.S.S. 

296.0 

21.3 

55040-60 

arXiv: 1103.0763 

PKS0447-439 

H.E.S.S. 

72.4 

-43.8 

55174-84 

arXiv: 1303.1628 
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6 Discussion 


Hadronic interactions predict neutrino emission in the TeV-PeV range associated with a flux 
of gamma rays. The prediction that the total neutrino energy flux F^, is approximately equal 
to the total high-energy photon flux is relatively robust, at least when attributing this 
emission to a 100% hadronic origin [48, 49]. As it is not affected by pair production losses, 
the neutrino emission is expected to have energies systematically higher than the electromag¬ 
netic component. The Fermi measurement of the source flux in the 0.1-100 GeV range is an 
underestimate of the overall electromagnetic spectrum, since it covers only three decades of 
energy. This underestimates, around a factor of 2-3 [50], has to be taken into account in 
the comparison between neutrinos and gamma-rays. Moreover, the lack of TeV observations 
makes a direct comparison at the >TeV range difficult, and necessitates an extrapolation 
of the LAT data over several orders of magnitude. This extrapolation is performed using a 
fit of the Fermi data with a log-parabola or a power-law function, both of which have been 
used by the Fermi Collaboration to build its source catalogues. To see how the neutrino 
upper limits and photon energetics compare, the gamma-ray spectral energy distributions 
(SEDs) of all sources have been produced using the SED builder of the ASI Science Data 
Centre (ASDC) [51] adding, if needed, VHE data taken from the literature, building a hybrid 
photon-neutrino SED [52]. The limits for different spectral indices, from -3 to -1, are extrapo¬ 
lated from the limits obtained in the E~'^ energy spectra using the corresponding acceptance 
curve. Figure 14 shows the value of the neutrino flux limits for the five brighest sources as a 
function of the spectral indices. The range in energy for each limit corresponds to the 5-95% 
range of the energy distribution for a given spectra. 

Figure 15 shows the hybrid SED for the lowest p-value source 3C279. The shaded yellow 
area represents an extrapolation of the flux during the studied flares from the average flux 
observed by Fermi. The lower bound is computed from the average flux during the 2008-2012 
period, while the upper bound is simply the renormalised flux according to the maximum flux 
measured in the light curve. With this simple criteria of the energy budget, the limit set by 
ANTARES for the blazar 3C279 is of the same order of magnitude as the electromagnetic flux 
measurements during the flares. It reinforces the need to search for a neutrino signal during 
the outburst periods when the gamma-ray flux and the accompanying neutrino flux are much 
higher. Therefore, with more data, ANTARES should be able to significantly constrain a 
100% hadronic origin of the high-energy gamma-ray emission. Fermi has reported some very 
intense outburst periods of 3C279 between mid 2013 and end of 2014 [55, 56], periods not 
considered in this paper. 

Figure 16 shows the hybrid SEDs for the four additional bright blazars. These sources, 
classified as FSRQs, have their gamma-ray flux suppressed at high energy. Only the brightest 
FSRQ objects are detected at TeV energies. On the contrary, BL Lac objects have high- 
energy spectral components, favoured for TeV gamma-ray detection, but are fainter than 
FSRQs. Both FSRQs [57, 58] and BL Lac objects [59] have been argued to be potential 
sources of neutrinos. This comparison between the gamma-ray flux and the neutrino flux 
limit provides an indication as to how to build an optimised source list for future searches 
with ANTARES and its successor KM3NeT [60]. 

7 Summary 

This paper discusses the extended time-dependent search for cosmic neutrinos using the data 
taken with the full ANTARES detector between 2008 and 2012. For variable sources, time- 
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dependent searches are significantly more sensitive than steady point-source searches thanks to 
the large reduction of the atmospheric background over short time scales. These searches have 
been applied to 41 very bright and variable Fermi LAT blazars and seven TeV flares reported 
by H.E.S.S, VERITAS or MAGIC telescopes. The most significant correlation was found with 
a GeV flare of the blazar 3C279 for which one neutrino event was detected in time/spatial 
coincidence with the gamma-ray emission. However, this event has a post-trial probability of 
67%, and is thus compatible with background fluctuations. Upper limits were obtained on 
the neutrino fluence for the selected sources and compared with the high-energy component 
of the spectral energy distributions computed with GeV-TeV gamma-ray observations. These 
comparisons show that for the brighter blazars, the neutrino flux limits are of the same order 
of magnitude as the high-energy gamma-ray fluxes. With additional data from ANTARES 
and with the order of magnitude sensitivity improvement expected from the next generation 
neutrino telescope, KM3NeT, the prospects for future searches for neutrino emission from 
flaring blazars are very promising. 
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Figure 9. Illustration of the flare-selection concept using the maximum-likelihood block algorithm 
with the blazar 3C273. (a) Gamma-ray light curve using the 2008-2012 Fermi data with thershold of 
100 MeV. (b) Evolution of the likelihood L (red dots) as a function of the number of added blocks 
(Nb)- The green dots show the evolution after including the prior V. The solid lines indicate the 
selected maximum, (c) Histogram of the flux convolved with a Gaussian distribution of the error 
on the flux. A Gaussian fit, shown by the black line, has been performed using data over the range 
(illustrated by the green line) between 60% (on the left-hand side) and 98% (on the right-hand side) of 
the maximum value. The value of the baseline is indicated by the red line, (d) Example of a selected 
flare of 3G273. Dots, the green and the yellow histograms display the raw data, the denoised light 
curve and the selected flare. The red dashed line corresponds to the 5cr threshold. 
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Figure 10. Skymap in galactic coordinates showing the position of the 41 selected Fermi blazars 
(red circles, see Section 4) and the seven TeV blazars (green circles, see Section 5) together with the 
ANTARES visibility (dark blue is maximal). 
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Figure 11. Results for the blazar 3C279. (a) Event map around the direction of 3C279 indicated 
by the green cross. The full red (hollow blue) dots indicate the events (not) in time coincidence with 
the selected flares. The size of the circle around the dots is proportional to the estimated angular 
uncertainty for each event, (b) Distribution of the energy estimator dE/dX in a ±10° declination 
band around the source direction. The red line displays the value of the event in coincidence with the 
flare in a 1° cone around the source direction, (c) Gamma-ray light curve (red dots) with the error 
bars (black) measured by the LAT instrument above 100 MeV. The green and purple histograms show 
the denoised light curve and the selected flare periods respectively. The red and blue lines display the 
times of the ANTARES events associated with the source during a flaring state and other events in 
a 3° box around the source position, respectively. 
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Figure 12. Upper limits on the neutrino fluence for the 41 studied Fermi blazars in the case of E~'^ 
(green), E~‘^ exp{—E/10 TeV) (blue), E~'^ exp{—E/l TeV) (yellow) and E~^ (red) neutrino energy 
spectra. The number in parenthesis after the name of the source in the x-axis indicates the total 
effective flare duration At during the studied period. 
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Figure 13. Upper-limits on the neutrino fluence for the seven studied TeV blazars in the case of E~'^ 
(green), E~‘^ exp{—E/10 TeV) (blue), E~'^ exp{—E/l TeV) (yellow) and E~^ (red) neutrino energy 
spectra. The number in parenthesis after the name of the source in the x-axis indicates the total 
effective flare duration At during the studied period. 



Figure 14. Extrapolation of the neutrino upper limit energy fluxes (integrated in the 5-95% energy 
range) from to other spectral indices for the five brightest Fermi blazars. 
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Figure 15. Gamma-ray SED of 3C279, observed by Fermi/LAT and MAGIC in January 2009 (red 
dots), January 2011 (green dots) and June 2011 (blue dots) [53, 54]. The yellow/green dots are the 
average flux with 2008-2010 data (2FGL). The observed data points are corrected for absorption by 
the extragalactic background light assuming z = 0.536. The dashed lines represents the fits performed 
on Fermi data using a log-parabola (LP) function. The shaded yellow area represents an extrapolation 
of the flux during the studied flares from the average flux observed by Fermi. The lower bound is 
computed from the average flux during 2008-2012 period, while the upper bound is the maximum 
flux measured in the light curve (xl3.3 the average flux). Finally, the coloured solid lines indicates 
the neutrino upper limits for different spectral indices (from E~^ in red to E~^ in blue). 
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Figure 16. Hybrid SED of the four additional bright Fermi blazars PKS1510-089, 3C454.3, 4C21+35 
and PKS0426. Points and lines are as described in the caption to Figure 15. 
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